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ABSTRACT 

Context. The Kepler object KIC 12557548 b is peculiar. It exhibits transit-like features every 15.7 hours that vary in depth between 
0.2% and 1.2%. Rappaport et al. (2012) explain the observations in terms of a disintegrating, rocky planet that has a trailing cloud 
of dust created and constantly replenished by thermal surface erosion. The variability of the transit depth is then a consequence of 
changes in the cloud optical depth. 

Aims. We aim to validate the disintegrating-planet scenario by modeling the detailed shape of the observed light curve, and thereby 
constrain the cloud particle properties to better understand the nature of this intriguing object. 

Methods. We analysed the six publicly-available quarters of raw Kepler data, phase-folded the light curve and fitted it to a model for 
the trailing dust cloud. Constraints on the particle properties were investigated with a light- scattering code. 

Results. The light curve exhibits clear signatures of light scattering and absorption by dust, including a brightening in flux just before 
ingress correlated with the transit depth and explained by forward scattering, and an asymmetry in the transit light curve shape, which 
is easily reproduced by an exponentially decaying distribution of optically thin dust, with a typical grain size of 0.1 pm. 
Conclusions. Our quantitative analysis supports the hypothesis that the transit signal of KIC 12557548 b is due to a variable cloud of 
dust, most likely originating from a disintegrating object. 
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1. Introduction 



Rappaport et al. (2012 1 recently reported the discovery of a pe- 
culiar signal in the light curve of KIC 1255754^] A transit-like 
signal is seen every 15.7 hours, with a depth greatly varying 
from event to event, ranging from less than 0.2% to more than 
1.2%. The transit timing is consistent with a single period and 
shows no variations down to the 10~ 5 level. This suggests that 
the changes in transit depth cannot be produced by a fast preces- 
sion of the orbit, which would in turn require the gravitational 
perturbation of a second companion in the system. In addition, 
the phase-folded light curve shows no evidence for ellipsoidal 
light variations down to 5 parts in 10 5 . This sets an upper limit 
of about 3 Mj up on the mass of the orbiting object. The expla- 



nation proposed by Rappaport et al. (2012 1 entails the disinte 



gration of a super-Mercury, caused by the high surface tempera- 
tures. Material stripped away from the planet surface forms and 
stochastically replenishes a cloud of dust that obscures part of 
the stellar disk, causing the dimming of the stellar light to ran- 
domly vary from orbit to orbit. Only a planet with mass on the 
order of Mercury would allow dust particles to accelerate to the 
escape velocity and form a trailing cloud of sufficient size. 

To test the interpretation of Rappaport et al. ( 2012| l, a quan- 
titative understanding of the observed transit light curve is re- 
quired, in particular the asymmetry of the in-transit portion of 
the light curve and the excess of flux before ingress. We re- 
duced the complete set of publicly available Kepler observations, 
phase-folded the light curve and grouped the nansits based on 



their depth, as explained in Sect. 2] A model for the trailing dust, 
described in Sect. [3] is then fitted to the data to estimate the tran- 
sit parameters, the basic shape of the cloud and the average size 
of the dust particles. The results of our analysis and their inter- 
pretation are discussed in Sect. [4] 

2. Analysis of Kepler data 

The star KIC 12557548 was observed by Kepler in long- 
cadence mode, delivering one frame every 29.4 minutes. This 
results in only 3-4- in-transit points per orbit, and therefore only 
an averaged and smoothed version of the light-curve can be re- 
covered by phase-folding the photometric time series. The pub- 
lic release of January 2012 extended the available data span to 
~ 1 .4 years of observations (six quarters) and includes more than 
22,000 photometric points. To combine and normalize the en- 
tire series, we designed a dedicated data reduction approach that 
starts with the raw Kepler photometry, which is obtained with 
simple aperture photometry. The Kepler pipeline also produces 
decorrelated photometry, using an automated decorrelation pro- 
cedure. This, however, often delivers over-corrected light curves, 
in which the astrophysical signal is reduced or even cancelled to- 
gether with the instrumental effect^] Therefore the decorrelated 
photometry data are not suitable for our purposes. 

A visual inspection of the raw light curve identifies the fol- 
lowing features: 

1 . Jumps in flux from quarter to quarter, at the level of several 
percent; 



Henceforth denoted with KIC 1255 b 



http://keplergo.arc.nasa.gov/ContributedSoftwareKepcotrend.shtml 
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Fig. 1. The normalized light curve of KIC 1255b, phase-folded 
using a period of P — 15.6854 hours. A binned light curve (solid 
squares, with error bars) and the best-fit model (solid line) are 
overplotted. 

2. Jumps or fast variations in flux within a quarter, at the per- 
cent level or less; 

3. A periodic modulation of the stellar flux, probably due to 
spots on the rotating surface or other intrinsic stellar vari- 
ability; 

4. Transits-like features with a well-defined cadence but vari- 
able depth. 

The instrumental systematics can be efficiently modelled 
with a linear combination of the Cotrending Basis Vectors 
(CBVs), which track the state parameters of the spacecraft and 
the instrumentation in timd^ We determined the optimum num- 
ber of CBVs for each quarter by minimizing the variance of 
the out-of-transit data, while visually checking that the periodic 
modulation attributed to stellar activity is preserved as much as 
possible. For each set of state vectors, we determined their co- 
efficients via linear regression and divided the data by the linear 
model. In this way different quarters of data were combined and 
instrumental systematics were minimized while preserving the 
astrophysical signal. To remove the low-order stellar modula- 
tion, we discarded the data points belonging to the transits and 
averaged the remaining points in each 15.7-hour orbit. This op- 
eration was performed separately for each section of contiguous 
data (i.e. data without temporal gaps) to avoid flux discontinu- 
ities. We then spline-interpolated the averaged points in time and 
divided the light curve by this baseline stellar signal. 

2.1. Transit light curve 

We phase-folded the normalized light curve using a period of 



P = 15.6854 hours, as determined by Rappaport et al. ( 2012| l 



The resulting transit light curve is shown in Fig. [1J with the in- 
dividual data points in light grey. A relative photometric preci- 
sion of 8.2xl0~ 4 per-point is achieved, which is ~30% higher 
than expected from pure photon noise. When binned by 0.008 in 
phase (black squares in Fig. [TJ, the scatter in the data is reduced 
by a factor of ~13, i.e. the square root of the number of points per 
bin, suggesting that systematic noise is not significantly present 



3 see http://archive.stsci.edu/kepler/release_notes/release_notes 1 2/ 
DataRelease_12_2011113017.pdf 



Fig. 2. Distribution of the transit depths, as measured by fitting 
the Kepler photometry with our model light curve with all pa- 
rameters but the transit depth fixed to the best-fit values. Orbits 
that don't sample the in-transit portion of the curve are dis- 
carded. 



in the data at this level of photometric accuracy. The binned light 
curve exhibits the following characteristics: 

1. An asymmetry between ingress and egress, with the latter 
being significantly longer than the former; 

2. Absence of a secondary eclipse within the photometric pre- 
cision; 

3. A small increase in flux just before ingress. 

To assess the dependence of the shape of the transit light 
curve on transit depth, we divided the events into two ex- 
treme classes: shallow (depth between 0.2% and 0.5%) and deep 
(depth > 0.8%). These ranges were chosen to maximize the dif- 
ference in depth while maintaining adequate photometric preci- 
sion. 

We determined the transit depths by first fitting the binned 
light curve, as explained in Sect. [3] to determine the overall pa- 
rameters of the transits. Subsequently, each separate orbit was 
fitted with the same curve, keeping all parameters fixed except 
for the depth, which we determined by estimating a scaling fac- 
tor via a x 1 minimization. The resulting distribution of transit 
depths is shown in Fig. [2] It peaks around 0.6% and extends 
from zero depth to about 1.3%. The negative transit depths in 
the histogram are due to noise in the photometry. The absence of 
a clustering around zero depth suggests that almost every orbit 
observed by Kepler shows a transit-like feature, except for ~8 
days during quarter 2, in which the signature seems to be absent. 
This results in an upper limit of the planet radius, as determined 
in Sect.|4] 

The light curves for the deep and shallow transits, shown in 
Fig. [3] share the same features of the global light curve except 
for the flux excess before ingress, which is not evident in the 
curve for the shallow events. The deep and shallow light curves, 
binned by 0.008 in phase, have a relative photometric precision 
per-point of 1.7xl0~ 4 and 1.2xl0~ 4 respectively. The solid lines 
in Fig. [3] show the best-fit model light curves, which were ob- 
tained by separately fitting the free model parameters. 
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Fig. 3. Light curves of KIC 1255b, obtained by dividing the events into two groups based on their transit depths and binning the 
photometry by 0.008 in phase. The left panel shows the average light curve of the deep events (deeper than 0.8%), and the right 
panel shows the average light curve of the shallow events (depths between 0.2% and 0.5%). The best-fit models for the two light 
curves are shown as a solid line. 




Kepler's view 



Fig. 4. Graphical representation of some of the quantities used in 
our model. The stellar and orbital radius are to scale, while the 
position of the orbiting body is shown with a black dot. 

3. The dust model 

We created a simple dust model with a minimum number of free 
parameters to quantitatively model the transit-like light curve. 
The most relevant variables are shown in Fig. [3] We assumed 
a cloud of dust distributed on a circular orbit at the same ra- 
dial distance as the parent body (a - 0.013 AU). The latter is 
supposed to not contribute to the eclipse in a measurable way, 
as suggested by the absence of any transit-like signals during 8 
days of Kepler photometry. It is therefore also reasonable to as- 
sume that the vertical extent of the cloud is negligible compared 
to the stellar diameter. Hence our model is one-dimensional, and 
it is expressed in terms of 0, the angle between the observer, the 
center of the star and the orbiting body. This translates to an or- 



bital phase tf according to = Imp. The planet crosses a chord of 
length r c on the stellar disk, corresponding to an impact param- 
eter of b — [1 - (r c /2R t ) 2 ]^ 2 for a stellar radius R*. The crossed 
chord, as seen from the cloud orbit, has an angular extension of 
f c = arcsin(r c /2a), which is a free parameter of the model. 

The cloud of dust is assumed to be optically thin, and its 
extinction cross-section p is expressed in units of stellar area. It 
drops exponentially away from the planet and is parametrized 
by a multiplicative factor c e and by an exponential parameter A, 
according to 



p(A0) = -^e- 
xRt 



A(&8) 



= c e e 



-A(M) 



(1) 



where AO — (6 - ff) is the angular distance between the position 
of the planet and an arbitrary point along the orbit, measured 
clockwise. The total model light-curve in units of stellar flux as 
a function of orbital phase ip is 



- f 

6 k 



[1 - / e (0) + uw\ d0, 



(2) 



i.e. the stellar flux corrected for the extinction component I e and 
the scattering component I s , convolved with a boxcar function 
that simulates the 29.4-minute exposure time of Kepler, which is 
long compared to the intrinsic temporal variations. The exposure 
time At translates into an angular range of 5 = 2nAt/P, where P 
is the orbital period of the planet. 

The extinction component I e {0) is given by 



4(0) 



/-*2tt 



p{6 - ff) i(ff ',r c ) dff ', 



(3) 



which is the convolution between the extinction cross-section 
and the intensity of the crossed stellar disk i(0' , f c ). The latter is 
described by the limb-darkening law i(p) — 1 - u(l - p), where 
p is the cosine of the angle between the line of sight and the 
normal to the stellar surface, and u - 0.79 is the linear limb- 
darkening coefficient appropriate for a K4V star in the V-band 
( Claret|2000 1. Expressing p in terms of ff and f c gives 



Ke f ,r c ) = l 



1 



a [ 

rZ^ 



sin 2 (r c /2) 



sin 2 ff 



(4) 
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set to zero for r c /2 < 6' <2n - f c /2. 

The scattering component I s (ff) is given by 



h(0) 



■In 



p(6 - 0) p(8') d0 



(5) 



where p{ff) is the Henyey-Greenstein (H-G) phase function 
( |Henyey & Greenstein|1941| i 



p(0) 



4tt(1 - 2gcos6' + g 2 )' 1 



(6) 



set to zero for \n-9'\ < r c /2 to account for the secondary eclipse, 
and containing a single free parameter, the asymmetry parameter 
g with - 1 < g < 1 . 

In Eq. <|5j we assume that the star is a point source. Since 
the scattering is beamed with a cone angle <p for an asymmetry 
parameter g, and g cos <p (Lamers e t al.|1997) , the effects of 
the stellar disk become relevant only when (f> % 21° , the stel- 
lar diameter seen from the orbit of the cloud. This corresponds 
to g 'Z 0.89 and therefore our best-fitting values (see Tab. [TJ 
should not be significantly affected by the assumption that the 
star is a point source. Therefore, we do not include the inte- 
gration over the stellar disk, which significantly speeds up the 
model computation. Finally, the single-scattering albedo w, i.e. 
the ratio between the extinction cross-section and the scattering 
cross-section, is a free parameter. 

The model is fitted to the binned light curves with a Markov- 
Chain Monte Carlo (M CMC) with a Metropolis-Hastings algo- 



rithm (Hastings 1970 1. Sequences of 5 x 10 steps are gener- 



ated starting from different initial conditions, and merged after 
discarding the first 10% of data to account for the burn-in. The 
mixing between the chains is tested ( Germa n" & Rubin|1992| l and 
the + 1 cr errors of the parameters are computed by measuring the 
interval between 15.8% and 84.2% of the merged distribution. 
The best-fit values and their uncertainties are listed in Tab. [T] 
for the average light curve and for the deep and shallow transits, 
respectively. 



Table 1. Best-fit parameters and their Icr uncertainties, as de- 
rived from the MCMC analysis. From top to bottom: impact pa- 
rameter, time of mid-transit (in BJD-2,455,000), decay factor, 
total extinction cross-section (in units of stellar area), asymme- 
try parameter and single-scattering albedo. 



Parameter 


Average 


Deep 


Shallow 


b 
To 
A 

Ce 

8 
m 


0.63+0.03 
451.9449(3) 

5.1+0.2 
0.030+0.005 
0.874!°*" 
0.65!™ 


46 +(, U2 

451.9461(5) 

6.4+0.3 

036 +0(xn 
o.oju_ 0002 

74.^+0-044 
u./tj 064 

>0.94 


0.61+0.04 
451.9430(7) 

3.86+0.25 
0.012+0.001 
0.87 (fixed) 
0.65 (fixed) 



Note that a similar forward-scattering contribution is present in 
egress, but it is masked by the dominating absorption of the tail 
of dust. Finally, the absence of a secondary eclipse is explained 
by strongly reduced backward scattering. 

4.1. Dependence on depth 

A comparison between shallow and deep transits shows that the 
total extinction cross-section decreases with transit depth, and 
that its exponential drop is slower for shallower transits. The 
photometric precision of the shallow light curve is insufficient to 
detect the forward-scattering peak, therefore we fixed the asym- 
metry parameter and the albedo to g = 0.87 and m = 0.65, 
respectively, as resulting from the fit of the global light curve, 
which has the best photometric precision. Although t& and g for 
the deep light curve significantly differ from those for the av- 
erage light curve, we cannot determine whether or not this re- 
flects different physical properties of the particles. In fact, in our 
MCMC simulations the extinction cross-section and the asym- 
metry parameter are moderately correlated, due to the fact that 
the measured transit depth is the summed contribution of the ex- 
tinction and scattering components. It is impossible to solve this 
degeneracy with observations at a single wavelength. 



4. Results and discussion 

The model parameters describing the fit to the average light 
curve do not exhibit typical signs of bias due to the fact that 
we are averaging events of all depths: i) different MCMC chains 
mix well, indicating that the solution we find is independent of 
the initial parameter estimates, and ii) the model parameters are 
intermediate between those of the shallow and deep transit light 
curves, as expected if the global curve is a combination of the 
two. Therefore, we will start with discussing the model parame- 
ters obtained for the average light curve. 

The transits of KIC 1255b occur at a moderate impact pa- 
rameter of b — 0.63. The asymmetric shape of the curve, with a 
sharp ingress and a slow recovery of the flux in egress, is well 
explained by a trailing cloud with a density that decreases away 
from the parent body. The characteristic length of the tail, I /A, 
is about 2.5 times the diameter of the stellar disk as seen from 
the cloud orbit. 

The H-G asymmetry parameter g being close to +1 indi- 
cates particles that strongly favor forward scattering over back- 
ward scattering. This reproduces the observed flux excess just 
before ingress. The fact that this happens slightly earlier than the 
ingress is due to the slow Kepler cadence and the limb-darkened 
stellar profile. Without implementing these two in the model, it 
is impossible to reproduce the position of the scattering peak. 



4.2. Constraints on particle size 

To constrain the particle size in the cloud, we compared realis- 
tic phase functions of silicates with the best-fitting H-G function 
obtained from the MCMC analysis of the average light curve. 
Following Bohren & Huffman] ([1983), we calculated the opti- 
cal properties of particles given their material (Draine & Lee 
1984), shape and size distribution. We assumed spherical grains 
(i.e Mie scattering is the dominant process), as expected if the 
grains liquify before evaporating due to the high temperature. 
The particle size distribution is described by a power law with 
a -3.5 exponent, with upper and lower limits left as free param- 
eters. The resulting phase functions were compared to a range 
of H-G phase functions constructed by varying m and g within 
their Icr uncertainties. A distribution of particles with sizes be- 
tween 0.04 and 0.19 pm best agrees with the explored range of 
H-G functions. Although our results are far from conclusive, the 
analysis suggests that the grains composing the cloud have an 
average size of ~ 0.1 pm. 



4.3. The size of the parent body 

As explained in Sect. |2.1| during 8 days of Kepler observations 
no transit-like feature is visible. The phase-folded light curve, 
constructed using the points in that time interval, has an in-transit 
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photometric accuracy of 2.8 xl0~ 4 per point. This translates into 
a lcr upper limit of 1.15 R® for the radius of the parent body. 

5. Conclusions 

A simple, one-dimensional model of a trailing dust cloud can 
reproduce the observed light curve of KIC 1255 b in exquisite 
detail. Features such as the flux excess before ingress, the sharp 
drop in ingress, the pointed shape of the flux minimum and the 
long flux recovery after egress are all signs of a trailing cloud of 
dust that periodically occults the star. The most likely origin of 
this trailing dust is a planetary or sub-planetary body that evap- 
orates due to the intense stellar irradiation. If this is indeed the 
case, we are presented with the unique opportunity to probe the 
interior of an exoplanet. Although the geometry of the transit is 
well constrained by the currently available Kepler data, further 
observations with a faster cadence and multi-wavelength cover- 
age are required to not only constrain the size but also the com- 
position of the grains. 
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